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ABSTRACT 

The cosixLological bound objects were considered to form from the local maxima of 
cosmological density fluctuations; often assumed to be Gaussian random fields. In 
order to study the statistics of the objects with hierarchical merging, we propose the 
skeleton tree formalism, which can analytically distinguish the episodic merging and 
the continuous accretion in the mass growth processes. The distinction was not clear 
in extended Press-Schcchtcr (PS) formalism. The skeleton tree formalism is a natural 
extension of the peak theory which is an alternative formalism for the statistics of 
the bound objects. The fluctuation field smoothing with Gaussian filter produces the 
landscape with adding the extra-dimension of the filter resolution scale to the spatial 
coordinate of the original fluctuation. In the landscape, some smoothing peaks are 
nesting into the neighboring peaks at a type of critical points called sloping saddles 
appears, which can be interpreted as merging events of the objects in the context of 
the hierarchical structure formation. The topological properties of the landscape can 
be abstracted in skeleton trees, which consist of line process of the smoothing peaks 
and the point process of the sloping saddles. According to this abstract topological 
picture, in this paper, we present the concept and the basic results of the skeleton tree 
formalism to describe (1) the distinction between the accretion and the merger in the 
hierarchical structure formation from various initial random Gaussian fields; (2) the 
instantaneous number density of the sloping saddles which gives the instantaneous 
scale function of the objects with the destruction and reformation in the mergers; (3) 
the rates of the destruction, the reformation, and the relative accretion growth; (4) 
the self-consistency of the formalism for the statistics of the mass growth processes of 
the objects; (5) the mean growth history of the objects at the fixed mass. 
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1 INTRODUCTION 

Hierarchical clustering scenario, including the cold dark 
matter (CDM) model, may be the most established one 
for reconstructing various observational properties in the 
cosmological structure from the galaxies to the clusters of 
galaxies. Press & Schechter (1974) firstly proposed an ana- 
lytical formalism which derives the number density of bound 
virialized objects of the mass at any given epoch, with the 
assumption that the primordial density fluctuations is ran- 
dom Gaussian field. The mass function predicted by the PS 
theory shows reasonably the good agreement with N-body 
simulations even if it has more low mass objects (e.g. Lacey 
& Cole 1994). To reconstruct the observational properties 
in theoretical galaxy formation scenario, there are also ap- 
proaches which study the history of the mass growth for 
bound objects and the characteristic times (e.g. Lacey & 



Silk 1991; Kauffmann, White & Guiderdoni 1993; Cole et 
al. 1994). Most of them were based on the extended PS for- 
malism, which was proposed by Bower (1991) and Bond et 
al. (1991). It can derive the number density of objects of a 
certain mass at a given time subject to a larger object at a 
later time. Using the formalism, Lacey & Cole (1993; LC) 
calculated the "merger" rate. 

The PS formalism, however, has a limitation for de- 
scribing the history of the mass growth about the individual 
objects. The "merging process" described with the PS ap- 
proach in LC, cannot be interpreted as the same meaning of 
the merger in astronomical sense, in which the objects lose 
their identity. In the mass growth history for the astronom- 
ical objects, the continuous accretion onto a bound object 
without the loss of identity has different meanings from the 
mass accumulation with the loss of the identity in the ma- 
jor merger. The formalism based with the PS approach can- 
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not imply the distinction between "tiny" and "notable" cap- 
tures. For solving this problem, Manrique & Salvador-Sole 
(1995, 1996) proposed a formalism named "ConflUent Sys- 
tem of Peak trajectories" (CUSP) formalism as an extension 
of adaptive windowing by Appel and Jones (1990). They can 
be categorized into a type of "peak" theory (Doroshkevch 
1970; Adler 1981; Bardeen et al. 1986, hereafter BBKS), 
which can count the number of the peaks of the density 
related to the collapsing threshold, applying a low-pass fil- 
ter of the bound object scale to the fluctuation field. The 
CUSP formalism, unfortunately, needs the iterative calcula- 
tion to estimate the destruction and the reformation rates 
with mergers. This point becomes a disadvantage when we 
try to apply this scheme to the semi-analytical studies for 
the galaxy formation, including the mass accumulation his- 
tory of bound objects with Monte Carlo method as shown 
in LC. 

In order to give analytical description for the destruc- 
tion and the reformation rate in the merger, we developed a 
new approach called the skeleton tree formalism, using the 
topological characteristics in the smoothing of the random 
field with Gaussian filter. With the appearance of sloping 
saddles in the landscape of the smoothed field, we can pick 
up the merging events, and distinguish the merger and the 
accretion. The topological feature in the landscape with the 
sloping saddles can be extremely reduced into skeleton trees. 

In this paper, we will focus on the concept and the ba- 
sic description of the skeleton tree formalism. The outline of 
the paper is as follows. In Section 2, in order to distinguish 
the accretion and the merging in our context, we sketch the 
topological characteristics in the landscape of the smoothed 
field with critical points; the peaks and the sloping saddles, 
define the merger events with appearance of the sloping sad- 
dles, and reconstruct the growth history of the objects with 
skeleton tree picture. In Sections 3, 4, and 5, we formulate 
the mathematical description of the constraints, the proba- 
bility distribution functions, and the scale functions for the 
critical points. In Section 6, we shows the results of the evo- 
lution rates with the accretion and the merging obtained 
from the skeleton tree formalism with its consistency. Fi- 
nally, we present our conclusion in Section 7. We have rele- 
gated the details of the derivations to five appendices. 



This Fourier transform is represented with the window func- 
tion which is the Fourier transforms of the selection function; 



F(k;i?) = S('k)W('k;R) 



(3) 



In our interest cases, the window function works as a low- 
pass filter. 

We shall restrict ourselves to isotropic homogeneous 
Gaussian random fields with zero mean as descriptions of 
the initial fluctuations. For the field, the power spectrum is 
then only a function of fc = |k|; \5k\^ — |5(k)|. The fiuctua- 
tion spectrum filtered with the scale R is 



P{k-R) = \Sk\'W{k-Rf 



(4) 



In this paper, we take a normalized isotropic Gaussian filter: 
W(k;i?) = M/(fc;i?) = c-*=''^'/\ (5) 

In the next subsection, we will discuss the reason related to 
its unique property for the smoothing. 

In the linear theory of gravitational instability for the 
structure formation, the amplitude of the field in the over- 
density area firstly grows in proportion to D{t), where D{t) 
is the linear growth factor. According to BBKS, a bound ob- 
ject collapses from the area of a comoving scale R when the 
density of a peak in the fluctuation smoothed over the reso- 
lution scale R exceeds above a flxed threshold Scfi- Instead 
of viewing the peaks to be growing in density amplitude 
relative to the flxed threshold Sc-o, we can interpret that 
the threshold level 5c is decreasing as 5c-oD{to)D{t)^^ with 
fixing the initial fluctuation field F{r; R), where Sc-o was de- 
termined from the threshold at the present time io- In this 
paper, we take Einstein-de Sitter model: f2o = 1, Ao = 0, in 
which the relative threshold level Sc/Sc-.o ~ D{to)D{t)^^ — 
(1 -I- z), where z is the redshift at the time t. 

For standard initial fluctuations like CDM models, in 
general, the rms of the smoothed fleld < F{r; R)^ > is de- 
creasing as increasing R. For such a fluctuation, as decreas- 
ing the collapse threshold with time evolution, we can pick 
the collapse objects in the larger scale, which gives a reason- 
able sketch of hierarchical clustering picture. We will relate 
the filtering process and the hierarchical clustering descrip- 
tion in the next. 



2 IDENTIFICATION OF ACCRETION AND 
MERGING 

2.1 Hierarchical Evolution from Fluctuations and 
Filtering Process 

We will express the density fiuctuation field as the functions 
of the comoving spatial coordinate r and k ; 



5{v) = / d^ke*'' '■<5(k) . 



(1) 



Interested collapse objects of a comoving scale R can be 
identified as peaks greater than a threshold, whose fields 
are smoothed with a low-pass filter of the resolution scale 
R. The fluctuation, smoothed with the selection function 
S(v; R), can be expressed as 



F{v-R) 



Jd''roS{ro;R)5{r-ro) 
/rf3r„S(ro;i?) 



(2) 



2.2 Landscape of the Fluctuation Field in 
Position and Resolution Space 

Consider the random fleld in one dimensional (1-D) posi- 
tional coordinate of x is smoothed with a low-pass filter of 
a resolution scale R. This field is reproduced as a landscape 
which extends in two dimensional (2-D) extended space of 
{x,R). A smoothing peak with a low-pass filter make a ridge 
which is running along the direction of increasing R in the 
landscape. The threshold level of the collapse can be inter- 
preted as an ocean surface which makes a "shoreline" and 
some "lakes" and "islands" in the landscape. A "cape" on 
the shoreline of 5c, can be identified as a bound object at 
Sc. Then, we can count the number of the bound objects 
of the scale R to pick up the capes at the R in the land- 
scape. As the level of the ocean surface of 5c is decreasing 
with the evolution of the universe, the shoreline moves to 
the offing. It means that the bound objects grow their scale 
continuously. 
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However, if "islands" or "lakes" appear in the land- 
scape, they confuse the identification of "capes" as bound 
objects. As an example, let consider an island in the offing 
of a cape on a shoreline. In this case, we have also another 
offing side cape of the island; there are two capes which can 
be counted as the bound objects of "clouds". In these situa- 
tions, the cloud of the island cape contains the smaller cloud 
of the cape on the shore. These problems appear when we 
take general filters except the Gaussian filter. Bond et al. 
(1991) had shown the filter dependency of the landscape ( 
see their Figure 1 and 2 ), in which they take sharp fc-space 
filter and Gaussian filter. In the case of the sharp fc-space 
filter, the ridges cannot decrease monotonically as the res- 
olution scale R increases. On the other hand, in the case 
of the Gaussian filter, the ridge decreases monotonically as 
the resolution scale R increases. They represent about this 
feature as there is no "lake" of finite extent, and the ocean 
shoreline have no bays. Another property of the Gaussian fil- 
tering is that the variance is also monotonically decreasing 
as increasing R. In general, this picture is also valid in 3-D 
random field. With the help of this feature, we can distin- 
guish the accretion and the merge in the landscape produced 
from initial fluctuation fields. 



2.3 Definition of Monotonic Accretion and 
Merging in the Landscape 

In order to understand the monotonic evolution of the field 
smoothed with the Gaussian filter of the resolution scale R, 
we rewrite the derivative of the field as 

9F(r; R) 



dR 



= i?VV(r; R) 



(6) 



This is identical to a diffusion equation of the variables {R^, 
r). For aU critical points of V^F < (V^F > 0) like peaks 
(holes) , it guarantees the monotonically decreasing (increas- 
ing) as the scale increasing as 



OFir;R)^^^^aF^^^^ 



dR 



OR 



(7) 



This monotonicity of the peak smoothing also guarantees 
that the peak runs continuously on the ridge to the shore 
cape without the island in the landscape. It can be rea- 
sonable that the smoothing and the merging of peaks are 
defined as the continuous accretion growth and the merging 
event of bound objects. If we have islands in the landscape 
with other filters of Gaussian, however, we cannot distin- 
guish the accretion and the merging with the confusion for 
the scale identification of the related bound objects. For- 
tunately, we can exclude this problem as long as using the 
Gaussian filter which guarantees the absence of island in the 
landscape as shown above. This is the reason that we take 
the Gaussian filter in this paper. 

A ridge in the landscape, then, represents continuous 
accretion growth of a bound object. On the other hand, 
some ridges terminate on the slope of neighboring ridges. 
The vanishing point of the ridge on the slope of the neigh- 
boring ridge can be defined as a type of critical points. We 
call it as sloping saddle since it is a saddle point on the 
slope of the neighbor peak. The sloping saddle can represent 
the reasonable feature that a bound object loses the iden- 
tity loss in the merger, associated with a tree structure in 
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Figure 1. The schematic presentation for the abstraction steps 
from the smoothing fluctuation field to the skeleton tree picture. 
(a) Schematic representation of the landscape. The symbols of 
and means the peak and the sloping saddle, respectively. 
The peak trajectory with filter smoothing makes the ridge. A 
ridge terminates at a sloping saddle at {Re, 5c) which is associated 
with a ridge of the neighboring peak through the resolution scale 
interval from R^ to Rj . (b) The critical points and the ridges are 
picked up as the abstract tree of the fleld around a sloping saddle, 
(c) The tree is presented by the reordering in (/?, 5) space. The 
hatched area means the merger, (d) The graphical presentation 
of the smoothing peaks with the merger and the accretion in the 
skeleton tree formalism. 



which the branches of the ridges are nested at the junctions 
of the sloping saddles. Then, we can translate the topology 
in the landscape to the tree structure named with skeleton 
tree, which consists of the accretion branches and the merger 
junction picked up with the sloping saddles. 



2.4 Skeleton Tree Picture 

We will describe the abstraction from the landscape to the 
tree structure as presented in Fig. 1 schematically. In Fig. 1 
(a), the local structure of the landscape is represented with 
contour lines, and the different classes of interest critical 
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points are marked; the peak as and the sloping saddle as 
(g), and the ridges are represented as the dashed lines. In this 
example, the sloping saddle appears on a resolution scale Re 
at a threshold Sc- A ridge terminates at the sloping saddle 
iSi{Rc, 5c). It means that a bound object of the scale Re loses 
its identity at the threshold Sc- There is a ridge neighboring 
with the sloping saddle, on which the peak has the density 
of 5s at the same resolution scale of iSi{Rc,5c). The peak 
on the neighboring ridge is reforming around the sloping 
saddle. Then, we can consider that the reformation starts 
from {Re, Sa) of a peak neighboring with the sloping saddle. 
The merger with the destruction and the reformation starts 
from 6a and ends at Si with the scale range from Ra to Ri 
in the case of Fig. 1 (a). 

As shown in Fig. 1 (b) as an abstraction of Fig. 1 (a), 
then, we can introduce a tree structure associated with the 
sloping saddle iSiiRcSc) and three peaks; the neighboring 
peak Q{Rc,Sa) of iSiiRcSc) in the same resolution Re, the 
progenitor peak Q{Ra,5a) on the same ridge of (8>(7?c,<5c) 
at Sa of the density of Q{Rc,Sa), and the reformed peak 
Q{Ri, Si) on the same ridge of the neighboring Q{Rc, Sa). 

For describing the history of the hierarchical cluster- 
ing with the merging, it is convenient to reorder the tree 
structure along the threshold level instead with the filtering 
scale. Because the threshold level is identical to the time 
as monotonically mapped with Sc — Sc-fiD{t)~^ . Fig. 1 (c) 
schematically shows this reconstructed feature from the Fig. 
1 (b). The continuous accretion growth is represented as the 
dashed lines along the ridges, while the merger is represented 
as the hatched area in Fig. 1 (c). In this case, the merger 
occurs during the interval between Si and Sa- Even if we can 
follow the neighboring peak on the same ridge over Sc in 
the landscape, however, we have two peaks at a threshold 
during the merger after the reordering with the threshold, 
as shown in the hatched area in Fig 1. (c). It means that 
the reforming peak on the ridge neighboring with the slop- 
ing saddle also loses its identity with the merger. Then, the 
neighboring ridge is divided into two part of a disappearing 
peak > Sc and a reforming peak < Sc- 

The merger feature is simplified as a joint, which con- 
nects three lines of two disappearing peaks and the reform- 
ing peak at a point of the sloping saddle as shown in Fig. 
1 (d), where the lines are the abstraction of the ridges as 
the continuous accretion growth. In the Figure, the line pro- 
cesses of the accretion branches are connected with the point 
process of the merger at the {Re, Sc), which makes a skeleton 
tree. 

In the derivation of the skeleton tree, we have neglected 
the detail features of three peaks around ®{Rc,Sc) except 
counting the number of the associated lines, and approx- 
imated as the destruction and the reformation occur at 
{5c,Rc) instantaneously. This is reasonable since no other 
parameter sets except {Sc,Rc) can be defined without any 
extra parameter for modeling the merger event. The consis- 
tency of this simple picture can be proved with the result in 
the section 6. 

We can practically reproduce accretion and merging 
with the script of skeleton tree as shown in in Fig 1(d); the 
arrowed line is related to the mass growth with continuous 
accretion, and the joint of hatching triangle is a merger event 
with the destruction and the reformation with the mark ® 
of the sloping saddle. We note that the branching number 



at the merging points is always two in the progenitor side. 
This is natural since the merger rate of three or more multi 
peaks are negligible to that of double peaks. 



3 THE CONSTRAINTS OF CRITICAL POINTS 

3.1 Line and Point Processes in the Skeleton Tree 

As shown in the previous section, the skeleton tree formalism 
is described by a set of the line and point processes of the 
smoothing peak along the ridge and the sloping saddle in 
the extended space (r, R). The line process of the smoothing 
peaks in (r, R) is equal to the point process in the original 
spatial coordinate r. 

Their density fields of these point processes are de- 
scribed as sums of S functions : 

npk{v;R) = ^5(=*)(r-rj,fc.,), (8) 

i 

n,,(r,i?) = ^5'^''\r-raa,i)S{R-Raa,^), (9) 

i 

where the subscripts of pk and ss means peaks and sloping 
saddles, Raa is the resolution scale of the sloping saddles. 
Since UpkCpv and risad^rdR are the numbers in 3-D infinite 
volume of d^r and 4-D infinite volume of d^vdR, we call the 
former and the later as the spatial density and the instan- 
taneous spatial density, respectively. 

We can express the point processes entirely in terms of 
the field and its derivatives with the spatial coordinate r 
and the resolution scale R- In the neighborhood around a 
critical point, with its constraint of VF(r)|cr = , wc can 
expand the field in a Taylor series: 

F(r) ~ F(rer)+ ^jV^VFjerArAr 



1 



' V® VF|crArArAr , 
and its derivatives can be also expanded as 



+ 3!^' 



VF(r) ~ V®VF|crAr, 
VF(r) ~ V®VF(r)|cr + V( 



' VFIcrAr 



(10) 

(11) 
(12) 



where Ar = r — Vcr and the suffix cr means the value at 
the critical point. 

The critical points can be divided into non-degenerate 
one and degenerate one. The extrema like a peak 
and hole can be categorized into non- degenerate one. 
Provided the condition of the non-degenerate extrema 
det(V ® VF)|cr / , Eq. (11) can be rewritten to 



= (V®VFj,0"'VF(r) 



(13) 



Using the second derivatives of the field, the number density 
of the extrema can be represented as 



= 5<^>((V®VF)"'Vi^) 
= |det(V® VF)|<5<^^(VF) 



(14) 



In order to describe the point process for the extrema, thus, 
it is enough to take the terms to the order of the second 
derivatives. 

In the degenerate case of V (8 VF|cr = , however, we 
cannot describe the displacement vector only with the first 
and second derivatives of the field as the non-degenerate case 



© 1999 RAS, MNRAS 000, ??-?? 



Statistics of Merging Peaks of Random Gaussian Fluctuations: Skeleton Tree Formalism 5 



of Eq. (13). The sloping saddle is a kind of the degenerate 
critical points. 

Under the transformation of the principal axis in the 
spatial coordinate, the part of six components related to the 
second derivatives V(8> VF in the covariance matrix becomes 
diagonal. It means that the second derivatives have three 
eigenvalues as 



(Fii,F22,F33) = -cr2(Ai,A2,A3) ,Fa,fj =0(/3/a) 



(15) 



where (T2 is the rms of V (^ VF ( see the definition in Ap- 
pendix A ). All the eigenvalues are not null in the non- 
degenerate cases. On the other hand, the degenerate critical 
points has one null eigenvalue at least as A3 — 0, where 
we assumed Ai > A2 > A3 for convenience. This is the 
reason for the break of the non-degenerate condition as 
detV ® V-F|cr = in the degenerate case. In general, a 
sloping saddle has a neighboring peak at the degenerate di- 
rection. In this degenerate direction, we cannot take the in- 
version of Eq. (11) as x:i—X3^cr = i^ssl^r^-Fs since Fsslcr = 0. 
In order to describe the point process of the sloping sad- 
dles, we use the expansions of F33 at the degenerate direction 
and V'^'F for the rest 2-D non-degenerate components in a 
couple of Eq. (12) and Eq. (11) as 



F33 ^ y^ F:i:ii3Axi3 ~ F333Aa;3 , 



p=i 



7(2) 



(2)«T7(2)paJ2) 



V^^^F ~ \/^^> (g)\/^^>FAr'-''> , 



(16) 



(17) 



where the suffice of (2) means the 2-D non-degenerate space. 
The set of the expansions derives 



J2) ^(2) 



= (V<^^®V<^^F1,0"'V<^^F(: 



X3 - X3,ss = -^333-^; 



333 -f^ 33 



(18) 
(19) 



Furthermore, we should remember that the sloping sad- 
dle is defined in the extended space with the resolution scale 
R. The gradient \/F{r) can be expanded with the derivative 
of i?; 

VF(r) ~ ^^\ss{R " Rss) = V^VFUs RiR - Rss) ■ (20) 
In the degenerate direction, this gives 



R — Rss — I / F3aaR 



F3 



(21) 



We obtain the constraint for the instantaneous spatial 
density of the sloping saddles as 

uss = ^5(=''(r-r,,,,)5(i?-i?s.,») 

i 

= |det(V(2>®V(2'F)|5(2)(V(2'F) 

3 

X \F333\5{F33) l^F3„ci?jc5(F3) . (22) 

a=l 

3.2 Constraints of Peaks and Sloping Saddles 

Under the transformation to the diagonal principal axis, the 
peaks require the condition of Ai > A2 > A3 > which 
should be added to the above condition of the extrema: Eq. 



(14). With the additional condition of the peaks, the con- 
straint for the spatial distributions of the peaks is 

3 3 

C(F''"'|pcaks) = ai\\i\2\3\Y[S{F^)Y[e{Xi) , (23) 

i— 1 a — 1 

where 9{Xi) is the Heaviside function. 

The sloping saddles requires the additional condition 
of Ai > (i = 1,2) as excluding the merger of the hole. 
With the constraint of Eq. (22), the constraint for the in- 
stantaneous spatial distribution of the sloping saddle can be 
described as 

C(F' 'I s. saddles) = o-2(J3i?|AiA2'u;3(wi + W2 + ^3)! 
3 2 

x<5(A3)[]<5(F,)[]e(AO, (24) 

where we use asWa = F3aa- 



4 PROBABILITY DISTRIBUTION FUNCTION 

The joint probability distribution of n-dimensional random 
variables with multivariate Gaussian can be described as 

exn[ip(") T/-.(nxn)-lp(n)l 
p(F(n))^p(n) ^ f^^PU£ ^ ^ ^ ^ ^ dF'") . (25) 

(27r)"|detC(")|V'' 

The convariance matrix c'"^"' is defined as the expectation 
value of the direct product of the vector F'"' : 

^(n>in) _ pi(n) p,(n) 



(26) 



where the suffix (n) of the F'"' is the dimension of the pa- 
rameter space. 

In general, the conditional probability of the event A 
with the constraint event B is given by the Bayes for- 
mula P(A|B) = P(A,B)/P(B). The joint probability, then, 
can be expanded with the conditional probabilities. If the 
parameters in the vector are all Gaussian distributed, we 
can directly obtain the covariance matrix of the condi- 
tional probability. There is a general theorem which is ex- 
tremely useful when we calculate the above joint prob- 
ability from the conditional probabilities. For the count 
of the peaks, we need only the 10 dimension parameters 
p(io) ^ (^p^ ^p^ y ^ yp) Pqj. ^Yie count of the sloping sad- 
dles, we should extend the parameter space to the 20 di- 
mension F^^"' = (F, VF, V ® VF, V (g) V (g) VF). The 20 
dimension covariance matrix can be found with the explicit 
forms of the parameters in Appendix A. With the help of the 
theorem (see Appendix Bl), the divided conditional proba- 
bilities are described for the case of the peaks and the sloping 
saddles in Appendix B2 and B3, respectively. 

For a type of critical points in a n-dimensional param- 
eter space described with constraints, in general, the prob- 
ability weighted density is 

ncr(F("')dF("' = P(F("')C(F("' I cr. points)dF("' . (27) 

Then, the ensemble averaged density of the type of criti- 
cal points restricted with m-dimension parameters, can be 
obtained from the integration over (n-m) dimension param- 
eters; 



< ric 



^("i)\ 



>dF 



(m) 



dF 



{n — m) 



n,r{F'-">)dF 



p{m) 



(28) 
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The ensemble averaged density gives the scale function we 
interest. 



5 SCALE FUNCTIONS OF CRITICAL POINTS 

We wiU derive the differential number densities of all the 
peaks, the nesting peaks, the non-nesting peaks, and the in- 
stantaneous differential number density of the sloping sad- 
dles, which are presented as the scale functions and the in- 
stantaneous scale function of Npk{R,5)dR, Nnest{R,5)dR, 
Npk{R,5)dR, and 'Nss{R,5)dRd5, respectively. 

The scale functions for all the peaks, the non-nesting, 
and nesting peaks are briefly given in the appendix CI, C2, 
and C3, according to BBKS and the CUSP formalism. We 
had checked that the contribution of the nesting peaks is 
negligible as shown in Fig. 2 . Thus, we treat the scale func- 
tion of peaks as that of the non-nesting peaks hereafter. 

In the first subsection, then, we sketch the derivations 
of the scale functions of peaks. In the next subsection, we 
present the new calculation for the instantaneous scale func- 
tion for the sloping saddles. 



5.1 Scale Function of Peaks 

The scale function of Npk{R,S)dR can be calculated from 
the ensemble-averaged density of the peak Npkiy, R)di' =< 
Upkiv'jR) > dv, where v = 5/ao, according to Manrique & 
Salvador-Sole (1995). The scale function means the differen- 
tial number density per infinitesimal range of 7? at a fixed 
i^, not from that per infinitesimal range of i^ at a fixed R. 
The forward one and the last one are denoted by a roman 
capital and a caligraphic capital which is the same as the 
notation of BBKS, respectively. The spatial density for the 
peaks of the non-degenerate critical points with a certain 
filtering scale can be calculated by the same way of BBKS. 
In order to obtain the density of peaks with up-crossing 
a certain 5c in the range of the filtering scale between 
R — AR and R, we should pick up the critical points which 
are equal to or smaller than Sc at the filtering scale of R 
and becomes larger than 5c at the smaller filtering scale 
of R — AR. It means that we should count the peaks of 
5c + {dFph/dR)AR < Fpk < 5c on scale R . As shown in 
the discussion of the Gaussian filtering, the condition of the 
counting is expressed as 



5c + V^FRAR <F <5c 



(29) 



As obtained for the peaks in BBKS, in general, we can 
obtain the spatial number density of the peaks in infinite 
ranges dvdx of the density contrast v and the second deriva- 
tive of it with a certain filtering scale R; Npk(v,x; R)dvdx. 
We can transform this number density to the number satis- 
fying the condition of the our counting; 



Npk{R,5c) = 



lim — — / dx I dvMpk{jy,x;R) 
Afl^o AR 7„ 7_^^^ 

Npk{uc;R)<x>pk ( ^4^ \R. 



MR) 



l^cb 



ao[R) 



(30) 
(31) 



where < x >pk is the mean value of x defined from the 
distribution function as 



< X >pk = 



Hpk jy, -yvc 
Gpk{'y,yvc 



(32) 



In practical treatment, we can use the averaged mapping 
relation 



duc ~< x >pk 



02{R) 

MR) 



RdR . 



(33) 



Fig. 2 illustrates Npk{R,5c)RRl and N''^'\R,5c)RRl 
versus R/Ro with the present collapsing threshold condition 
5c — 1.68 in the power-law fluctuations of n = 0,-1, and 
—2. The former and the latter are represented with the solid 
and dotted lines, respectively. The critical scale of Ro, was 
determined by 5c/oo{Ro) = 1 with 5c = 1.68. The number 
of the nesting peaks is negligible compared with that of all 
peaks. Hereafter, then, we will take N{R,5c) = Npk{R,5c). 
The PS scale functions, obtained from the same filtered ran- 
dom fields, are also presented with the dashed lines. We also 
present the two scale functions of peaks for standard CDM 
(SCDM) model with the normalization of 8 Mpc with the 
bias factor b — 1. Even if the two functions in the small scale 
are different from each other between the peak theory and 
the PS formalism as pointed out by Appel & Jones (1991), 
they deviate relatively little over the cosmological interest 
interval ( 2 decades in the filter scale, 6 decades in the mass 
scale around Ro). 

Fig. 2 shows that the scale functions of the peaks 
Npk{R,5c)RRf) become proportional to R^'^ asymptotically 
in the small scale. As shown in BBKS, the cumulative num- 
ber is a useful quantity which can be evaluated analytically: 



npk{vc = — oo) 



l/c^ — OCi 



^'*"-iw«-'- <«' 



As i? ^ 0, the density contrast Vc ^> Q and the peaks of the 
low contrast Vc dominates in the cumulative number. Then, 
Npk{R,5c)RRo oc npk{uc — > 0;i?,) ex npk{vc = — c»;i?,) in 
the small scale. Since i?* is proportional to the resolution 
scale R, the asymptotic feature seen in Fig. 2 is not unex- 
pected thing. 



5.2 Instantaneous Scale Function of Sloping 
Saddles 

The instantaneous scale function of sloping saddles can be 
directly calculated as its ensemble-averaged density: 



^ss{R,Vc)dRdVc 



< ■n.ss{i'c\R) > dRdVc 



(35) 



The brief calculation is described in the appendix C. 

We should note that the scale function is defined as 
the density in the infinite interval dRdVc. The scale function 
of the peaks Npk{R,5c)dR introduced before is the differ- 
ential density of peaks per infinitesimal ranges of the reso- 
lution scale i? at a fixed density threshold 5c. On the other 
hand, the instantaneous scale function of the sloping saddles 
Nss(i?, 5c) is the differential density per infinitesimal ranges 
of R and 5. Their difference come from the fact that the 
point processes of peaks and sloping saddles can be defined 
in the 3-D spatial coordinate and the 4-D space extended 
with R, respectively. Then, the instantaneous scale function 
is responsible for the time evolution properties of the scale 
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Figure 2. The scale functions are presented for various fluctuation models of the initial density fluctuation. The solid and dotted lines 
represent the scale function of the all peaks and the nesting peaks. The dashed line is that of the PS formalism with the same fllter. 



function with the merging process as seen in the skeleton 
tree; the Une of the smoothing peak connected with other 
lines at the junction of the sloping saddle. In the next sec- 
tion, using the skeleton tree formalism, we will present the 
evolution of the scale functions with the distinction between 
the merger and the accretion according to the skeleton tree 
formalism. 



6 EVOLUTION WITH MERGER AND 
ACCRETION 

6.1 Instantaneous Scale Functions for the 
Disappearing and Reforming Peaks 

As discussed in the sketch of the skeleton tree picture with 
Fig. 1, the reformation of a peak can be only identified with 
the destruction of two peaks around a sloping saddle. In 
the above section, we can obtain the instantaneous number 
density of the merger of a scale R at 5c in the infinite interval 
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of d6c- Thus, with changing between the variables (5c and Vc, 
the instantaneous scale functions of the destruction and the 
reformation for the peaks of the scale R at 5c are described 
as 



N'^{R,Sc)dRdSc 
N^{R,Sc)dRdSc 



2TSiss{R,5c)dRdi^c -, 
Nss{R,Sc)dRdPc, 



(36) 

(37) 



where the superscripts of / and d mean the reformation and 
the destruction in the merger, respectively. 

With the formula, the net change of the differential 
number density of the peaks in dR during the interval of 
dSc is simplified to 

S{R,5c)dRdSc = N'^{R,5c)dRd5c-N'^{R,5c)dRdSc 

= -Nss{R,5c)dRduc. (38) 

It means the net destruction of the peaks is the same as the 
appearance of the sloping saddles. This is reasonable since 
the formation of the peaks does not occur in exact sense and 
all peaks only disappear at the sloping saddle in the original 
filtering process. From the concept of the skeleton tree, the 
peak neighboring with a sloping saddle is interpreted to lose 
its identity with the merger. We considered that the ridge of 
the neighboring peak is divided into two part of a disappear- 
ing peak > Sc and a reforming peak < 5c, even if the ridge of 
the neighboring peak continues over Sc- Thus, the reforming 
number of the peaks can compensate for the destruction of 
them, and the appearing number of the sloping saddles lefts 
as the net disappearing number of the peaks. The consis- 
tency of this formula can be checked with the conservation 
equation for the scale function of the peaks, as shown in the 
next subsection. 



6.2 Instantaneous Rates for Accretion, 
Destruction and Reformation 

We will consider the instantaneous scale growth rate for a 
peak of R as the continuous scale growth of the accretion. 
The accretion growth rate can be described with the scaled 
Laplacian x as 



Racc{R,X,Sc)dt 



xdiR 



-dSc 



(39) 



We should note that the growth rate for a peak depends on 
its particular value x. According to Manrique & Salvador- 
Sole (1996), the mean growth rate for the objects of the scale 
R can be expressed as 



Racc(R,t)dt 



1 



xa2R 



> dSc 



1 



< X >pk (J2-R 



dSc , (40) 



where < > means the average of the function with x and we 
used the relation 



< X 



l^dxx-^ Npk{R,5c,x)dR 

J^dxNpk{R,Sc,x)dR 

^^dxNpk{vc,x;R)^JR 

J^ dx x J^pk {vc, X, R)^dR 

Gpk{-f,'y'^c 
Hpk{'y,^Vc 



--<x> 



pk 



(41) 



From this mean scale growth rate, we can define the relative 
growth rate with the accretion as 



r"-{R,t)dt 



dRaccjRX) 

dR 



(42) 



With the similar consideration applied for the merging 
objects, we can introduce the mean scale growth rate and 
relative growth rate contributed with continues accretion 
during the merger phase; 



^acc,my-^^ i^jyJ'T^ 



r^{R,t)dt = 



1 



< X >S3 O2R 

dRacc,m\R, t 



d5c 



(43) 
(44) 



where < x >ss^ Hss{'y,n,'yvc)/Gss{'y,K,'yi^c) ■ 

The instantaneous destruction and reformation rates 
can be defined directly from their instantaneous scale func- 
tions as 

N'^{R,S)dR 



r{R,t) dt 



N{R,t)dR 



dS 



(45) 



Thus, the conservation equation for the scale function 
N{R, t) is given as 



dN{R,t) ^ d{RaccN{R,t)) ^ ^^^ ^^ 



dt 



dR 



(47) 



where S{R,t) can be given as the net source term with the 
destruction and reformation rates of r'^{R,t) and r-' {R,t) ; 

S{R,t)^ [rf{R,t)-r''{R,t)]N{R,t) . (48) 

The conservation equation can be rewritten to 



d\nN{R,t) 
dt 



^ r^ {R,t) - r'^{R,t) - r''{R,t) , 



(49) 



which can be rewritten to 



^}^l^^l^=r%R,t)+rf{R,t)-r''{R,t)-r''{R,t),{50) 



where 



r=(i?,t) = -i?a 



din N{R,t) 
'' dR 



(51) 



is the rate of shift in the scale space with the accretion for 
the number density distribution. 

The inverses of these rates of r""'^ •'''''' {R, t) give the time 
scales of the individual processes. In Fig. 3, the rates for 
the reformation, the destruction, the relative growth with 
accretion, the relative growth with continuous accretion for 
merging objects and the shift in the scale with the accretion 
at the present are presented with solid, dotted, short dashed, 
long dashed, and long dash-dotted lines. The left hand side 
and the right hand side of the conservation equation (52) are 
also presented with short dash-dotted line and crosses, the 
former d In N{R, t)dt can be calculated with the partial time 
derivative of Npk {R, t) and the latter is calculated from the 
set of the individual rates as r''{R,t) + r^{R,t) ~r''-{R,t) — 
r-{R,t). 

We can see that the partial time derivative of 
dlnN{R,t)/dt can be reproduced well from the right hand 
side of the conservation equation from the individual calcu- 
lations of r%R,t),r^R,t),r''{R,t), and r°(R,t) as shown 
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Figure 3. The rates for the evolution processes at the present with Sc = 1.68. The soUd, dotted, short dashed, long dashed, and 
long dash-dotted lines mean the reformation rate r^\dt/dSc\, the destruction rate r'*|c(i/d<5c|, the relative accretion rate r°'\dt/dSc\, the 
relative accretion rate r'^\dt/d5c\ during the merger, and the shift rate r'ldt/dSd, respectively. The evolution rate of the scale function 
\dlnNpk{R,t)/dtdt/d5c\ and the right hand side of the conservation equation (51) is presented by short dash-dotted line and crosses. 
The absolute values arc presented for the two last rates. 



in Fig. 3. Especially, around the scale of 7? ~ Ro, both of 
them shows the same feature, in which the values shapely 
switch from negative to positive. Then, we can verify the 
consistency of our formalism with the equality between the 



left hand side and the right hand side of the conservation 
equation. 

The rates of the reformation and the destruction have 
a maximum around i? ~ 7?o and rapidly decrease for the 
larger scale. On the other hand, the accretion dominates 
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Figure 4. The evolution of growth rates at fixed masses for the 
SCDM model in an Einstein-de Sitter universe. The solid and 
dashed lines present the merger and the accretion growth rate 
t. In the figure, the highest curves at the right 



rf ■ t and 
is for M/Mo 



10" 



, 10^\10 



10 

-0.5 



, and successive curves are for M/Mq 
1,10°'^ and 10^ 



the merger processes of the reformation and the destruction 
and the rate becomes almost constant for the large objects 
( i? >> i?o )• It is consistent with the picture that the large 
object are growing in the cosmic time scale. These mean 
that the larger objects of R > Ro grows with the accretion 
in the linear regime of the gravitational instability, while the 
smaller objects are accumulated into the larger objects with 
the merging in the non-linear gravitational growth regime 
under the critical scale of i?o, which is related to the present 
threshold 5c. In the theory of gravitational instability, the 
smaller than i?o can become in non-linear growth. We should 
remark that the merger process is the most efficient around 
the critical scale i?o even if the merger still dominates the 
accretion in the small scale. These properties of the dom- 
inating merger around Ro are already suggested from the 
previous N-body calculations (e.g. Navaro, Frenk, & White 
1997 ). 

From the view point of the mass growth history for a ob- 
ject of the scale R, the results in Fig. 3 should indicate that 
the dominant growth process switches from the accretion to 
the merger accumulation around the threshold oi 5 = ao{R) 
and the merger accumulation rate becomes maximum and 
decreases rapidly. We can see directly this feature in Fig. 4, 
which represents the evolution of the growth rate with the 
reformation and the accretion at fixed scales. From these 
results, in general, the growth of a halo firstly starts with 
the accretion process, secondly switches to the merging and 



is suppressed soon after the merger dominates. This feature 
also was remarked from the simulations. 

We should remark that the results in Fig. 4 are not the 
mean growth histories for the individual objects since the 
fixed scales cannot be directly related to the final scales of 
the objects at the present due to the successive destruction 
and the reformation cooperated with the accretion. In or- 
der to reconstruct the mean growth history for a individual 
objects, we should extend the basic of the skeleton tree for- 
malism including the background effects. This extension of 
the skeleton tree formalism will be described in the following 
papers. 



7 CONCLUSION 

We have derived an analytical expression for the statistics 
and the evolution of the cosmic bound objects with the refor- 
mation, the destruction and the accretion. It is applicable to 
any hierarchical clustering models in which structure grows 
via gravitational instability. The skeleton formalism is de- 
rived as a natural expansion of the peaks theory of BBKS. 
In the landscape reproduced from the random field with the 
Gaussian filter, we have followed the smoothing of the peaks 
as the accretion growth of the objects, and have picked up 
the "sloping saddles" as the merging events with the destruc- 
tion and the reformation of the objects. The line and point 
processes of the peaks and sloping saddles in the landscape 
produce the tree structure. Then, we call our scheme as the 
skeleton tree formalism. With the skeleton tree formalism, 
we can estimate the rates of the reformation, destruction, 
and the accretion in any hierarchical clustering models. The 
set of these rates can reproduce the evolution of the scale 
function of the objects with the conservation equation. This 
reproduction of the evolving scale function verifies the self- 
consistency of the skeleton tree formalism. With the rate cal- 
culation of the individual processes, we can find the merger 
processes are efficient around the critical scale of Rq deter- 
mined as (5c = ao(-Ro)- The dominant growth process of the 
objects switches from the accretion to the merger accumu- 
lation around the critical threshold related to the scale. It 
is important to reproduce the mass growth history with dis- 
tinguish between the accretion and the merger when we try 
to reproduce the cosmic structure and the galaxy formation 
in the hierarchical scenario. 
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APPENDIX A: THE COVARIANCE MATRIX 

We will introduce 20-dimension vector F*^"^ = {F, VF, V ® 
VF, V (g) V (g) VF, ). The VF, V ® VF, and V (g) V (g) VF 
have three, six and ten independent components, respec- 
tively. They can be expressed in terms of Fourier transforms: 



F(r; R) 
F4r;R) 

Fc.fi(v;R) 

Fap^{r;R) 



d^ke*"'F(fc; i?) , 



dF{v;R) 

dXa 

d'^Fjr; R) 

dXadXjj 



d^e'^^'ko^Fik^R) 



d^ke*"'fc„fc^F(fc;i?) 



d'^F{r;R) 
dxadxfjdx-y 



d' ke' '^kakpkyF{k;R) 



(Al) 
(A2) 

(A3) 

(A4) 



It is useful to introduce the integrals over the filtered fluc- 
tuation spectrum: 



a;(i?) =47r 



dkk 



2j+2 



P{k-R) 



(A5) 



when transforming the above values to the non-dimensional 
ones. 

The covariance matrix C' ' is defined as the expecta- 
tion value of the direct product of the vector F' 



,(20) 



C. 



(20) 



.< F<'°>F 



(20) 



> 



We can represent the matrix as 
Mja 

(^(20) 

Mn 
M22 

M22 
M33 





/ 1 1/3 1/3 

1/3 1 1/3 
1/3 1/3 1 




4 



/M3 










(A6) 

(A7) 

(A8) 
(A9) 

(AlO) 
(All) 



M33 = 

M02 = 

Mi3 = 

^i _ 

13 — 



■ 1 1/5 1/5 ^ 

1/5 1/5 1/15 

,1/5 1/15 1/5, 



-? 



-i 



/j^(lx3) q(1x3)\ 



(Mi3,Mi3,Mi3,0' 



(3xl)>, 



(A13) 
(A14) 



M 



T 

6; ( 



(1, 1/3, 1/3) , e, = {5u, 52^,53^) , (A15) 

where k("^'"',I("^"' and O^"^™' is the (n x m) dimension 
matrix with every entry unity, the (n x n) dimension unit 
matrix and the (n x m) dimension null matrix, respectively. 



APPENDIX B: CALCULATIONS OF JOINT 
PROBABILITIES 

Bl Theorem for the Joint ProbabiHty Expansion 

According to Adler (1981) and BBKS, all Gaussian dis- 
tributed parameters represented by the n-dimension vector 
Z — (Y, X), which can be devised into m-dimension vector 
Y and (n-m) dimension vector X, imply us that the condi- 
tional probability F(Y | X) is a Gaussian 



P(Y 1 X) 



exp[-iAY < AY (g) AY I AX >-^ AY] 



(Bl) 



(27r)™/2| dot < AY ® AY I AX > |i/2 
where AY = AY- < AY | AX >, with the mean; 

< AY I AX > = < AY ® AX >< AX^ >"^ AX (B2) 
<AX^>"' = <AX®AX>"\ (B3) 

and the conditional covariance matrix: 

C(YIX) =< AY(g) AY I AX > 

=< AY (g) AY > 
- < AYg)AX>< AX^ >"^< AX® AY> , (B4) 

where < > represents the mean value and AX = X— < 
X >. With the theorem, we can obtain the conditional prob- 
abilities from the covariance matrix. 

B2 Divided Conditional Probability 

We can expand the joint probability distribution with the 
conditional probabilities as 



P{F, F^,F„fi, F^fi^)dFdF^J>dF';>dF^ 



(3)wp(6)^p{10) 



F{F^p^\F,F^,F^p)dFl^^] 

X P{F,F^,F^f,)dFdF^^^^dFf^^ 
P{F,F^,F^p)dFdFi^'^dF^^^ = 

P{F^0\F, F^)dF^^^P{F, F^)dFdF^''> 



(B5) 



(B6) 



Using the theorem above, these conditional probability func- 
tions are represented in the explicit forms: 



PiF,F^) 
P{F^p\F,Fc.) 



3^/^exp[-|Q(F,F„)] 

P{Fo.p\F) 

exp[-lQ(F^^|F)] 
(27r)3|defC(F„3iF)|i/2 



(B7) 



(B8) 
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P{F^fj-,\F,F^ 


,fa 


fi) = P(F„^^|F„) 










exp[-iQ(F,^^|F, 
(27r)5|deiC(F,/3^|F„ 


)|l/2 


, (B9) 


Q{F,F^) 


= 


9 9 ' 




(BIO) 


Q{Fo.p\F) 


= 


fI^C{F^^\F)-^ F^p , 




(Bll) 


Q{Fc.m\Fc.) 


= 


Fal3-fC{Fai3-y\Fa) Fa0^ , 




(B12) 


Fc^fi 


= 


F^f3 - Mo2a,7'F , 




(B13) 


-t' af3j 


= 


Fc,f3-y - Sa^^MnF^ , 




(B14) 



and the conditional covariance matrix 

C(F<,;3|F) = M22 - MiTaatT'Mo: 

^2 

M202) , 



^(M22 



M202 

C{F^^-y\F^) 



5at 



j^(3x3) q(3x3) 
q(3x3) q(3x3) 



M3 



M?;3ai"'Mi3 



(B15) 
(B16) 
(B17) 



B3 Probability Function for Peaks 



For the peaks, the original covariance matrix for the prob- 
ability is represented in the 10 dimension vector F^^'^'. As 
presented by BBKS, we can see that the three degrees of 
freedom related to the directional dependence drop off from 
the parameter space for the homogeneous fields. After intro- 
duce the eigen coordinate, we can represent the peak number 
density in 7 dimension parameter space. 

For convenience, we will introduce a new set of variables 

(B18) 
(B19) 
(B20) 



aov 


= 


F, 




CTlI/c 


= 


F,., 




0-23; 


= 


-vV = 


= Ai + A2 + A3 



(72^ = 



dF 



(J2Z 

afi 



(Ai - A3) 

2 
(Ai - 2A2 + A3) 

2 
|(Ai-A2)(A2-A3)(Ai 
X d vol[S'0(3)] 



(B21) 
(B22) 



A3)|dAidA2dA3 



= \2y{y^ - z^)\-aldxdydz d vol[5'0(3)] , (B23) 

where d vol[S'0(3)] is the volume element of the three- 
dimensional rotational group 5*0(3). We had used the re- 
lations : 



Ai 
A2 
A3 



-j{x + iy + z) 



(72 , 
(72 



2z) 



\{d\r 



(x - 3y + z) 



—a2dxdydz 



(B24) 
(B25) 
(B26) 

(B27) 



The probability function in these variables of 
{F,FcFap) is represented as that in the variables of 
{u,rj,x,y,z) with 77 = {1^1,^2,1^3)- 

P{F,F^,F^fj)dFd''F^d''Fo,fj = 



P{v, f), X, y, z, a, (3, 'y)dvd rj 
X \2y{y^ - z'^)\-a2dxdydz d vol[50(3)] 



(B28) 



where q, 13, 7 are the Euler's angles. The directional depen- 
dence is not so important in the isotropic field. Then the 
mean of peak number density is enough to consider the 
statistics of the density field. The mean can be obtained 
by the angle independently integrations. 



d vol[S'0(3)]P(i/, r), X, y,z,a, (5, 7) 



2n' 



27r^ 
= -^P{x,y,z\v)P{i^,rj) 



(B29) 



where we used J d vol[SO(3)] = 2-k^ /3\ with the degenerate 
of the triad orientation of the axis for the eigenvalues. 

The determinant and the inverse of the conditional co- 
variance matrix are given explicitly as 



|detC(..,|.)|^/^ = ^^^^-^(l-7^) 



2n1/2 



(72(70 



where 



1—7^ 

-151(^^3) 

/ 1 -1/2 -l/2\ 

R = 5 -1/2 1 -1/2 . 
V-1/2 -1/2 1 / 



(B30) 

(B31) 

(B32) 



The probability functions P{x,y, z\i'),P{iy) are the 
same of that described by BBKS as 



P{u,rt) = 

P{x,y,z\i^) = 

Q{x,y,z) = 



3^/' .1.2 



2 exp[--(!/ +377-77)] 



(27r)2 -'-''' 2 
33 . 55/2 



(27r)3-2- (1-72)1/2 

X exp[--Q(x,77,z)] , 

{x-x-tf , .,(. 2 , r 2 

■ + 15j/ +02 ,x. 



(1-7^) 



(B33) 

(B34) 

(B35) 



B4 Probability Function for Sloping Saddles 

For the sloping saddles, in general, we should start from the 
original covariance matrix in the 20 dimension vector with 10 
variables of the third derivatives of the field added to that for 
the peaks. After introduce the eigen coordinate, however, we 
need the third derivatives in only degenerate direction as the 
additional variables; v^aa = F^aa/us. As similar as shown 
in the case of the peak, finally, the statistics of the sloping 
saddles of the homogeneous field can be determined in the 
parameter space {i^,ri,x,y,z,'w), where w — {101,^12, wg) = 
(2^311, f322, t'333). The conditional probability with the third 
derivatives of w can be represented as 



P{W\U3) 
Q(w|!/3) 



exp[ 



■hQi^W^)] 



(27r)3/2|c«eiC(wJK3)!i/2 

W C(w|l^3)W , 



-3 w 

-—KM.ViU3 

5 



(B36) 

(B37) 
(B38) 
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53/2 . 3 . r. 



A 



(B39) 



(73(71 

The determinant and the inverse of the conditional covari- 
ance matrix for the variables V^cia are represented as 



detC{w\v'A)\^^'^ 


- 53/2.3-7^^ ""^ 


' 


(B40) 


C(wji.3)"' 


^ / Cll Co 


Co \ 
Co 
C33/ 


, (B41) 


(1 - ^2) -" C2. 
^ ' \Co Co 


Cll 


= 10 - 7«:^ , 




(B42) 


C22 
r.„ 


= C33 = 45 - 42k^ , 
21^2 - 15 




(B43) 
rB44^ 



With the condition of the first derivatives Fz for the 
sloping saddles, we can represent the exponent part of the 
conditional probability: 






Q\(w\\w2,W-},,Vz) = 



W\ — 



3(15- 14k^) 

(1-^2) 

(5 - 7k2) 



Q2(W2 1^3,1^3) 



Q-i{'Wz\v'>) = 



2(15 - 14k2 

3-5-7(25-21k^) 
22(15-14^2) 

_ (5 - 7^2) _ ■ 

™^- (25- 21.2) "^3 



W3 + TO2 



(B46) 



r'^aaa 



/i 



(25-21k2) 
We will list up some integrations as 

dwi exp — -Qi(wij«;2, W3,i^3) 
2^(1 -K^)^/^ 

31/2(15- 14.2)1/2 ' 

1 



I2 = 



^22 



dw2 exp 



(52(W2|W3, V'i) 



2V2^(15 - 14k 



2>,l/2 



(3- 5- 7)1/2(25 -21.2)1/2 ' 
1 



(B49) 



(B50) 



dw2 W2 exp 



2(w2\w-i,V-i] 



(5-7.^) _ 
dw-iY^ ~-, — ^^3 exp 



(25-21.2) 



--Q2(w2|w3,i^3) 



(5-7.^) 

(25-21.2) 



W3/2 



(B51) 



As shown in the definition of the sloping saddles, we are in- 
terested in mainly the probability for wz in the third deriva- 
tive components. The conditional probability function of W3 
can be obtained from the integration with the other third 
derivative components w\ , W2 as 



P{wz\v'i)dwz = / dw\ / dw'2,P(yf\vz) dwz 
© 1999 RAS, MNRAS 000, ??-?? 



2(27r)»/2(l-.2)i/ 


-xh 


X/2 




X exp 


-^Qim) 


dwz 




5 • 7I/2 






^27r(25 - 21.2) 




X exp 


-^Qim) 


dw-i , 


(B52) 


(25 - 21.2) (^^ + 


|..3) = 





(B53) 



Q{wz) = 



We can see the reasonable normalization in the integration 
of P{w3\vz) with Wz as 



dwzP{wz\uz) 



(B54) 



For the convenience, we will represent the probability of W3 
with vz = 0; 



(B45) P{mWz = 0) = 



—= — exp 

(25-21.2) 
52~7 



i™| 
2al 



(B55) 



(B56) 



We will calculate the integration required in the calcu- 
lation of the instantaneous scale function of sloping saddles; 





/ d' w\wz 


(toi + W2 + W3) -P(W 1^3) = 








<wi >ss (l + qss) , 


(B57) 


(B47) 




|.oo 






< wl >ss = 


/ dwziwzfPiwzli^z) 




(B48) 


= 


J-00 

(25-21.2) 
52.7 ' 


(B58) 




Qss = 


(5-7.2) 
(25-21.2) ' 


(B59) 



where we used 

dwiwi I dw2-P(w|i^3) dwz 



dw2W2 / dw\P{-w\vz) dwz ~ 

qssWzP{wz\uz)dwz . (B60) 



APPENDIX C: CALCULATIONS OF SCALE 
FUNCTIONS 

CI Ensemble Averaged Density of Peaks 

The constraint of the peaks can be rewritten to 

3 
C(F("''ipcaks) = \\i\2\z\Y[5{Fo.)e{\z) 

Q, = l 

{x-2z)[{x + zf-y']^a2,z 



33 



i-r 

(71 



'[[S{iyo,)eix-2z) . 



(CI) 
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With the constraint, in order to derive the ensemble- 
averaged number density of peaks, we will start from its 
probability weighted number density: 

npk{i^,x,y, z; R)dvdxdydz 

= 6 • P{v, Tj — 0,x, y, z)C{F ' \]ieaks)dvdxdydz 

5"^'3'/^ , 02 ,3 '^M-hA exp[-iQ(:r, y, z)\ 



(27r)3 cri ^1-72 

X '4)pk{x,y,z)(l)pk{x,y,z)dudxdydz , (C2) 

33 _, 

tljpk{x,y,z) = — (72 ^AiA2A3(Ai - A2)(A2 - A3)(Ai - A3) 

= {x-2z)[ix + zf-{3yf]yiy'-z') (C3) 
1 if f > y > and 

-y> z>y 

(t>pk{x,y,z) = { 1 iff>J/>fand 

3jy — a; > 2 > 2/ 
otherwise 

where 4>{x,y, z) is the condition of Ai > and we had mul- 
tiplied the probability expression by 6 to account for the 
ordering of the eigenvalues of Ai . 

After the integration over y,z, the ensemble-averaged 
density of the peaks is given as 

exp[— i;^^] 
< npk{v,x;R) > dvdx = dv 

(III ) Hi 

. fpkixf-^^^t(i^^£ll^z:Ql,, , (C4) 

where _R, = \/3— and 

CT2 
32^5/2 /•a:/4 



fpk{x) = 



dye ^^^ ''^ / il)pk{x,y,z)dze ^^ ^^ 



('x/2 py 

+ / dye-i^'''/^ / i>pk{x,y,z)dze-^''/^] 



x/4 

{x-^ - 3x) 



3y-x 



2 {Erf[(5/2)^/=^^] + Erf [(5/2)1/^ |]} 

^ l/2,,3]£^ 8 -5,2/g 

^^57r^ '^ 4 ^5^ 

+ C^-l)e-'^''']- (C5) 

The density of the peaks of 1/ at the scale R can be 
represented as 

exp(-ii/2) 



< npk{v\R) > dv = ^^^^/p3 ' Gpk{-y,Xt)du , (C6) 

where 



{2tiYRI 
Gpk{-1,Xt)= I dx fpk{x) , diy .[C7) 

Jo VMi-i") 

The function Gpk{'y, x,) have a following fitting formula ob- 
tained by BBKS, 

It;" - 37^^; + [Bw^ + Ci] exp[-^w^] 
C'pfcd.a:^*) - r~T^ r~Fi — i i^^) 



1 + C2 exp[— Csw] 



A = 
B = 



5/2 



(9-572) ' 
432 



(107r)i/2(9_ 57,2)5/2 



(C9) 
(CIO) 



Ci = 1.84 + 1.13(1-7^)^-^' 
C2 = 8.914- 1.27 exp(6.5l7^ 
Cs = 2.58exp(1.057^) . 

We also have the averaged x ; 



(CU) 
(C12) 
(C13) 



Hpk(ci,x,) 



X dx fpk (x) 



^ exp[-(x-.:,)V2(l-7^)] 
\ v/2^(l-7^) J 

= <x>pkGpk{-y,x,) . (C14) 

We can easily see that the mean value < x >pk is 
Hpk{'Y,Xf)/Gpk{'y,x»), which can be also represented as a 
approximated form as derived by BBKS : 



< X >pk = 



Oi + 02 cxp [-7/2(7iy/2)2] 
[6^ + 0.45 + (71^/2)2] 1/2 + ^^12 ' 
3(1 - 7^) , 6)2 = (1.216 - 0.97^) . (C17) 



(C15) 
(C16) 



C2 Scale Functions of Peaks 

The calculated ensemble-averaged densities is equal to them 
of BBKS as 



< npk{v, x; R) > dxdu — Npk{v, x; R)dxdv . 



(CIS) 



With transformation of the density contrast v to the reso- 
lution scale R, we can obtain the scale function of peaks: 

Npk{R, X, 5c)dRdx — Mpk{v, x; R)x ( — ) RdRdx 

{2tvYRI ao{R) 
^ exp[-(^^_^£/2£-7!)l^^.^ 

VMi-i^) 

This is the scale function of the peaks with the parameter of 
X. With the integration of x, the density function of peaks 
satisfying the condition of collapse threshold Sc is given as 

Npk{R,S,)dR = f xdxJ^pk{u,,x,R)^^RdR 
Jo ^o{R) 

Hpk{-f,X*) 1 2iCr2(-R) „ ,„ 

— exp[— -i^cl — TTTzRdR 



{2-kYRI 
= Mpki^cR) < X >pk 



2 "Vo(i?) 

02{R) 



oo{R) 



RdR. (C20) 



C3 Scale Functions of The Nesting and 
Non-Nesting Peaks 

The above peak counting does not distinct nesting peaks and 
non-nesting peaks. In order to count the nesting peaks, we 
will introduce the joint scale function of the nesting peaks : 

Npk (Rb, Rs', ScjdRtdRs ~ Npk {Rs,\Rb',Sc)dRs 

X f(Rb,Sc)dRt(P2l) 

where the fraction of the volume occupation with the objects 
related to the filtering scale Ri is given as 



f{R,5)dR = Vol{R) N{R,S)dR^ 



(C22) 
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M{R) 



Yol{R) = 



(C23) 



where N{R, S)dR is the scale function of non-nesting peaks. 
In general, the conditional density of the peaks at 5s 
within the background at 5t on the larger scale Rb can be 
expressed as 



N;r{Rs,5s\Rb,5b)dRs = 
Afp^{Rs,5s\Rb,5b)^4^ <^>Rs dRs 



(C24) 



in terms of the conditional density J\fpk{Rs,5s\Rb,5b) ob- 
tained by BBKS. Here, the function < x >=< x >ph (7,2^*) 
, where 



7^^ , 7 = 7 



7(1: 



^ 1 - e2 ' 

Vs(— j - (^Vb 

1 — ri 



(C25) 
(C26) 



The subscripts b, s and h mean the values of the scales Rb, 
Rs, and Rh {Rb > Rh > Rs )■ According to BBKS, we 
introduced the following notations; 



ri 



r2 = 





e - 


o'abO'os 


(C27) 


< k^ >h 

<k^>s 


^2 ^2 


O'2hO'0a 
O'0hO'2s 


(C28) 


<k^>h 

<k3>s 




O'shO'Os 
O'OhO'Ss 


(C29) 



ajh = 4:TT I dkk 



L2i+2 



5k rWik;Rb)Wik;Rs 



(C30) 



Since we use the Gaussian filter, we can rewrite as aj h ~ 

aj{Rh) with the rms average Rh — \{Rb + -Rs)/2] 

This conditional density of peaks gives a form for 



iV, 



pk 



{Ra,\Rb;Sc)dRs at the the limit of 5b 



In this paper, we use this simple form for the conditional 
density as 



iV, 



pk 



'{Rs,\Rb;5c)dRs^ 



AfpkiRs, 5c\Rb, 5c)^454 < ^ > -R. dRs 



(C31) 



The number density of the non-nesting peaks N{R, 5c) 
can be obtained by subtracting the density of the pairs of 
nesting peaks N"^"'^{R,5c) from that of all peaks Npk{R,5c), 
which includes the nesting peaks N"'^''^{R,5c) at 5c on the 
scale R. Then, the notations to the corrected scale function 
is 



N{R,5c)dR 



Npk{R,5c)dR~N'"''\R,5c)dR. (C32) 



The scale function of the nesting peaks can be calcu- 
lated from the integration of the joint scale function of the 
nesting peaks N^-'^^R, R'; 5c)dRcdSc, 



\R,5c)dR = / dR'N"^'\R',R;5c)dRc 
Jr 

dR'f{R',5c) 



N" 



xN;^'\R\R';5c)dR. (C33) 
This is a Volterra type integral equation, which can be solved 



with the iteration process in exact mean as shown in Man- 
rique & Salvador-Sole (1995). The number density of the 
non-nesting peaks N{R, 5c) can be also obtained with the 
same iteration process. The number density of the nesting 
peaks can be negligible to that of all peaks, as shown in Fig. 
2. 



C4 Ensemble Averaged Density of Sloping 
Saddles 

The constraint of the sloping saddles are rewritten to 

C(F<^"'|s. saddles)di? = 

1 Cr20-3 , 



\(x + 3y + z)(x — 2z)wi(wi + W2 + W3) 



3 cr2 

X 5(z ~ 3y + x)5'^'-''> (ri) ( 



RdR 



(C34) 



With the constraint, the probability weighted number den- 
sity of the sloping saddles is represented as 

iissiv, X, y, z, w; R)dvdxdydzd wdR 
= 6 • P{w\u;i = 0)d wP{x,y, z\v)dxdydzP{i',r] — Q)dv 
X C(F<^"^ls. saddles)dii 

o5/2 , c5/2 

-{^^)Ri'ss{x, y, Z)(j)ss{x,y, z) 



(27r)3 
exp[ 



1,,2 



X 



v'^]e:^p[-\Q{x,y,z)] 



dudxdydz 



X |w3(wi -f «;2 -I- 'u;3)jP(w|i^3 = 0)d w 
X (—\ RdR , 

Tpss{x,y,z) = y{y^ ~ z'^){x + 3y + z) 

X {x - 2z)5{z -3y + x) 
(j)ss{x,y,z) = (l)pk{x,y,z) , 



(C35) 



(C36) 
(C37) 



where the factor 6 is needed to account for the ordering of 
Xi the same as that of the peaks. With the integration over 
y, z, we can obtain the density of the sloping saddles similar 
to that of the peaks; 



< nss{i^,x,w; R) > dvdxd wdR 
exp[ 



-^u^]du 



{27T)2Ri 

^2^(1- 7=^) 

X \W3{wi +W2+ Wa)\P(w\u3 = 0)d^W 

—) RdR , (C38) 

where we used the relation {(J20-3/af)^''^ = \^/k^'^ /R,, 
and 



Jss\f^,X) — yiss\f^) 



x/2 

K) / dyy{y;x) 

x/4, 

/8r 



A..Me-^-/«[^(l- 



-15a;^/8x 



1 2,3x' 1 



(C39) 
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Ass(k) 



34 . 55/2 



'J 3y — x 

= 18y\x-2yf{4y-x) 

The density averaged with w is given 

exp[— i(!/^)]dz^ 
{2-kYRI 



(C40) 



(C41) 



< \iss(v,x;R) > dRdxdv 

X fss{x) 






da; 



X < u;2 >,, (1 + g,,)£4^Mi? . 

(Tl(R) 

With the integration over x, the averaged density is 
< nss{iy;R) > dRdv 



(C42) 



(27r)2E» 



dv 



(T3{R) 



X <wi, >ss (1 + qss)-^RdR , (C43) 



where 



fss(xy 



exp[ 



"2(1-7")- 



Tm^ 



dx . (C44) 



■7 



The exact result of the integration is obtained as an analyt- 
ical form 



Gss{n,i<.,xt.) = 



Ass (^J 



v/27r(l-7^) 



2(1-7^)2:., 12a;2 (39-5572), fi^ 

5(9-572)8^9-572)^ 5 ^ 



/27r(l-72) 24xt 4(27 - 3572)(9 - 572)r 



(9-572)9' 5 ' 25 

(783 - I23O72 + 5757'')(9 - 572)2 



+ 



750 



} 



X {1 + Erf 



375 V (6-572) ^ 



(9-572)(l-72) 

Xt 



jg 2(9-57^) 
} 



^2(6-572)(l-72) 



X e 2(6-57")] . (C45) 

The mean value of x at the sloping saddles is 

Hss{-y,n,x.t) 



<X>s 



Gss{'y,ii,x*) 



Hssi'y,K,x,) 



exp[ 



(x — x^) ■ 
"2(1-7")- 



x /ss (a;) ■ , 

^2^(1- 72) 



(C46) 



dx , (C47) 



where we can also obtain the analytical result of the inte- 
gration : 



Hssi'yjK.jXt) = 



Ass [f^J 



v/2^(l-72) 



4(1-72) 24:xt 2(99-11572) 



4{ 



+ ■ 



(9 _ 5^2)3 L 5(9 _ 5^2)2 25(9-572) 

9(1-7'^' -• 



6 — 572 



-}£ 



2 / 27r(l-72) 48x:; 8(57 - 6572)2^; 

^ 5 V (9 -572)7 "^(9 -572)2 + 5(9 _ 5^2) 

(2403 - 495O72 + 2675'r^)x, - 



75 



-} 



X {1+Erf 
2 /27r(l-72) 



375 V (6 - 572)3 



V (9-572)(l-72) 
{H-Erf 



}x, e 2(9-57") 
} 



V2(6-572)(l-72) 



X a:, e 2(6-57") ] (C48) 



APPENDIX D: SOME FORMULA IN THE 
CASES FOR POWER LAW SPECTRUM 

When the power spectrum can be represented as a simple 
power law as 

I 5fe |'= Ak" , (Dl) 

we can obtain the useful relations for the calculations as 



1 + 3 



a,{R) = ^i/2V2^i?-(^+:')r(!l±l ^ .V/2 



+ J) 



it* ^ 

2 

7 = 



(2i^LZi±l)!]_ij-2(«-.) 



2*-J(2j + n + l)!! 
6 .1/2 



ri 



r2 = 



'-R. 



n + 5 

_ai_ 2 ^ {n + i) 

0-002 (n + 5) 

ol .2 .n + 5, 



)2 = (::_L^) 

RsRb_-.(n+3)/2 



K 

Rs x2 
i?J 
^a n3 
i?J 



(D2) 
(D3) 
(D4) 
(D5) 

(D6) 

(D7) 
(D8) 
(D9) 
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